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ABSTRACT 

We present early-time photometric and spectroscopic observations of super- 
nova (SN) 2009kr in NGC 1832. We find that its properties to date support 
its classification as Type Il-linear (SN II-L), a relatively rare subclass of core- 
collapse supernovae (SNe). We have also identified a candidate for the SN pro- 
genitor star through comparison of pre-explosion, archival images taken with 
WFPC2 onboard the Hubble Space Telescope with SN images obtained using 
adaptive optics (AO) plus NIRC2 on the 10-m Keck-II telescope. Although the 
host galaxy's substantial distance (~26 Mpc) results in large uncertainties in the 
relative astrometry, we find that if this candidate is indeed the progenitor, it is 
a highly luminous {M v = — 7.8 mag) yellow supergiant with initial mass ~ 18-24 
M Q . This would be the first time that a SN II-L progenitor has been directly 
identified. Its mass may be a bridge between the upper initial mass limit for 
the more common Type Il-plateau SNe (SNe II-P) and the inferred initial mass 
estimate for one Type Il-narrow SN (SN Iln). 

Subject headings: galaxies: individual (NGC 1832) — stars: evolution — super- 
novae: general — supernovae: individual (SN 2009kr) 
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Introduction 



It is not yet exactly clear how to map massive stars of a given mass range to a core- 
collapse supernova (CC SN) subtype. We now know with a growing degree of confidence that 
solitary stars in the range of ~ 8-1 6 M inevitably explode as Type II-plateau supernovae 
(SNe II-P; e.g..lSmartt et al. 1120091 ). with much of their hydrogen envelope still intact (see 
Filippenko 1119971 for a discussion of SN classification). We also have evidence that the Type 
II-narrow (Iln) SN 2005gl had a luminous (My = —10.3 mag), very massive (possibly with 
initial mass M- m \ > 50 progeni tor star that exploded while in the luminous blue variable 



phase ( iGal-Yam fc Leonard II2009I ). This progenitor, the only example so far identified for 
a SN Iln, likely had a far smaller fraction of its outer H layer remaining than do SNe II-P. 
Various indirect clues may indicate that at least some of the Type Ib/c SNe are connected 
to the Wolf-Ray et phase, which is expected to occur for stars with M ini > 25-30 M (e.g., 
Crowther 1 I2007I ). That leaves the Type II-linear (II-L) SNe with no directly known pro genitor 
star, as well as the range M ini = 18-3OM (exactly the range that lSmartt et al. 1120091 dubbed 
"the red supergiant problem") without a well-established endpoint. From the rates derived 
with the Lick Observatory Supernova Search (LOSS) for a sample of well-studied SNe (Li et 
al. 2010, in prep.), we know that the majority of massive stars end their lives as SNe II-P, 
while the incidences of SNe II-L and SNe Iln are comparatively rare (~ 7% and ~ 6% of all 
CC SNe, respectively). 

Here we exam i ne the case of SN 2009kr in NGC 1832. SN 2009kr was discovered by 
Nakano &: Itagaki I (120091 ) on 2009 Nov. 6.73 (UT dates are us ed throughout) and was 



spectroscopic ally classified as a S N Iln (|Tendulkar et al. 1120091 ). and then as a "young 



type-II" SN (jSteele et al. 1 120091 ) . iLi et al. I (120091 ) first identified a possible progenitor 



star in archival HST images from 2004, using as reference a combined 160-s r' image from 
CFHT+MegaCarrl^l on Nov. 21.49. Here we show that both early-time photometric and 



1 Based in part on observations made with the NASA/ESA Hubble Space Telescope (HST), obtained 
from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., under NASA contract NAS 05-26555; the 6.5-m Mag- 
ellan Clay Telescope located at Las Campanas Observatory, Chile; various telescopes at Lick Observatory; 
the 1.3-m PAIRITEL on Mt. Hopkins; the SMARTS Consortium 1.3-m telescope located at Cerro Tololo 
Inter- American Observatory (CTIO), Chile; the 3.6-m Canada- France-Hawaii Telescope (CFHT), which is 
operated by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the 
Centre National de la Recherche Scientihquc of France, and the University of Hawaii; and the W. M. Keck 
Observatory, which is operated as a scientific partnership among the California Institute of Technology, the 
University of California, and NASA, with generous financial support from the W. M. Keck Foundation. 

10 A joint project of CFHT and CEA/DAPNIA. 
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spectroscopic observations strongly suggest that this object is a SN II-L, and we increase our 
confidence in the progenitor identification via further comparison of the HST images with 
Keck-II/NIRC2 adaptive-optics (AO) data. 



2. The Early-Time Nature of SN 2009kr 



2.1. Photometry 

Optical BVRI images of SN 2009kr were obt ained with the Lick Ob servatory 0.76- 



m Katzman Automatic Imaging Telescope (KAIT; iFilippenko et al. Il200ll ) and the 1.0-m 



Nickel telescope, and SMARTS +ANDICAM at CTIO They were all initially reduced using 
standard proce dures (see, e.g., iPoznanski et al. 1 120091 ) . Because of our follow-up campaign 



on SN 2004gq flModiaz 1 120071 ). also in NGC 1832, we had a well-calibrated photometric 
sequence in the host-galaxy field and te mplate images in all passband s for image subtraction. 
We used an image-reduction pipeline (jGaneshalingam et al. 1 120101 ) to reduce all data and 
calibrate the photometry to the standard Johnson BV and Cousins RI system. 

Near-infrared (NIR) observations were obtained with the 1.3-m Peters Automated In- 
frared Imaging Telescope (PAIRITEL; iBloom et al. 1120061 ) and reduced using standard pro- 
cedures. We used archival 2MASS images of NGC 1832 as templates to subtract from the 
SN images using HOTPANTS^] and calibrated the JHK S photometry against 2MASS stars 
in the field. 

The BVRIJHKg light curves are shown in Figure Ufa), relative to 5-band maximum 
(15.95 ± 0.02 mag on Nov. 13 ± 1, or JD 2,455,149 ±1). For comparison, we also show 



the light curves of the SNe II-L 1979C (Balins kaia et al. 



Barbon et al. |[l982bj ), 1980K flButa 1 11982 ; iDwek et al. 



980; de Vaucouleurs et al. 



1981 



1983 ). 1990 K (ICappellaro et al. 



19951). and 20~01cy ( unpublished KAIT data); the SNe I I-P 1999em JHainuv et al. 



2001 



(IPoznanski et al. 



Leonard et al. UgOOgJ) an d 1992H flClocchiatti et al. 111999 ): and,the possible SN II-L 2000dc 



20091 ). SN 2009kr does not fo llow the SN II-P pla teau; instead, it more 



closely resembles the SN II-L linear decline (see iBarbon et al. II19791 ). A steeper decline is 
present in VRI after ~ 65 days in the SN 2009kr light curves; however, this decline was also 
observed for SN 1980K in V. From our adopted distance and total ext inction fsee 32. 2D . 
SN 2 009kr reached M(5 max ) = -16.48 ± 0.30 mag, typical of SNe II-L jYoung & Branch 



19891 ). We thus conclude that SN 2009kr displays the photometric behavior of a SN II-L. 



http://www.astro.washington.edu/users/becker/hotpants.html. 
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2.2. Spectroscopy 



Low-resolution spectra of SN 2009kr were obtained on 2009 Nov. 10 an d 25, Dec. 9 and 



18, a nd 2010 Jan. 20 with the Lick 3-m Shane tel escope+Kast spectrogr aph ( IMiller Stone 



19931 ): on Jan. 8 with the Clay telescope+MagE (Marshall et al. 1120081 ): on Feb. 15 with the 



10-m Keck-II telescope+DEIMOS; and on Mar. 9 with the 10-m Keck-I telescope+LRIS 



The spectra were reduced and calibrated using standard procedures (e.g., iMatheson et al. 



20001 ) . The observing conditions were not photometric, resulting in un certain absolute-f lux 
calibrations. We observed with the slit placed at the parallactic angle (IFilippenko Ill982l ). 



We show the rest-frame sp ectral sequence in Fig ure [TYb); for comparison, we also display 



spect ra of the SN II-L 1980K (IBarbon et al. Ill982al): the SNe II-P 1 999em ([Leonard et al. 



20021 ). 1992H (IClocchiatti et al. Ill996) and 2004et flSahu et al. 1120061 ): and the SN II 2001cy 



at similar epochs. iPoznanski et al. I ( 120091 ) rejected SN 2001cy from their SN II-P sample, 
implying that it likely may be a SN II-L. Indeed, the SN 2001cy light curve is very similar (but 
of inferior quality and coverage) to that of SN 2000dc shown in Figu re [TTa). The followin g 



Leonard et al. 



Maguire et al. 



2002), 



20101) 



corrections for E(B — V)tnt have been made: SN 1999em (0.10 mag 
SN 1992H (0.09 mag Iciocchiatti et al. Ill996h . SN 2004et (0.41 mag 

SN 1980K (0.40 mag. lBarbon et al. Il982ah . and SN 2009kr (0.08 mag; see bel ow). SN 200 ley 
has b een corrected only for Galactic reddening, E(B — V)ghi = 0.21 mag (ISchlegel et al. 
1998h . 



The spectral sequence for SN 2009kr is typical of many SNe II: the earliest spectra 
exhibit a blue continuum, with relatively weak spectral features, followed by the onset of 
increased Balmer-line emission, together with the emergence of P-Cygni-like features. How- 
ever, clearly SN 2009kr differs from the canonical SN 1999em or the peculiar SNe 1992H 
and 2004et; in particular, the SN 2009kr P-Cygni Ha profile is dominated by the broad 
emission component. Several narrow emission lines appear superposed on the spectra, but 
these most likely originate from a neighboring H II region, seen ~ 1-2" northeast of the SN 
in an Ha image of the host galaxy (unpublished ESO archival data; the region's ionizing 
cluster can also be seen in the HST images). The SN 2009kr spectra more closely resemble 
those of SN 1980K or SN 200 Icy, which also exhibit relatively weak Ha absorption. The 
optical photometric decline rates [> 1 mag (100 day)" 1 ] of SNe 1980K, 2001cy, and 2009kr 
suggest that they are all SNe II-L (see § 12.11) . Since spectroscopi c features of SNe II appear 



to be correlated with their photometric behavior at early times (ISchlegel Ill996l ; iFilippenko 
19971 ). the similarities in the spectra of these three SNe are not unexpected. 



We also obtained a high- resolution optical spectrum of SN 2009kr on 2009 Nov. 25.46 
with CFHT+ESPaDOnS, with 4 exposures of 833 s each, from which we confirm that the nar- 
row Ha emission likely originates from the H II region. Moreover, we also used this spectrum 
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to estimate the reddening toward SN 2009kr via measurement of the Na I D equivalent width 
(EW) at the host-galaxy redshift (z = 0.006). We found EW(Nal Dl A5890) = 0.044 ±0.003 
A and EW(Nal D2 A5896) = 0.032 ± 0.004 A. Using the relation b etween extinction and 
EW(Nal D) from Elias-Rosa et al. (2010, in prep.), and assuming the lCardelli et al. I (119891 ) 
reddening law with updated wavelengths and a Galactic foreground E(B — V) = 0.07 mag 
(Schlegel et al. 1998), we derive E(B - V) tot = 0.08 ± 0.01 mag (E[V - J] tot = 0.11 ± 0.01 
mag), which we adopt for the SN. This relatively low extinction is consistent with both the 
SN color comparison (§ 12. ip and the overall blue continua seen in the early-time SN 2009kr 
spectra (Fig. [TJb]). 



3. Identification of the Progenitor Candidate 

Pairs of HST images of NGC 1832 were obtained with WFPC2 in bands F555W (~ V; 
460 s total) and F814W (~ J; 700 s total) on 2008 Jan. 11 (program GO-10877, PI: W. Li), 
as SN 2004gq follow-up observations. The SN 2009kr site is located on the WF3 chip (CX'l 
pixel -1 ). Cosmic-ray hits were rejected, and a 1600 x 1600 pixel mosaic of all four WFPC2 
chips was constr ucted using the STSDAS package routines crrej and wmosaic within IRAFlE 



Li et al. I (120091 ) were able to isolate the SN location in the WFPC2 images to 0.43 pixel 



(0"043) through comparison with a ground-based, post-explosion CFHT+MegaCam image. 

We were then able to better confirm the identification of this candidate through K p - 
band NIRC2 "wide" camera (0'.'04 pixel" 1 , 40" x 40" field of view) images obtained on 2009 
Nov. 28 with Keck-II+AO. Each of the 10-s frames was sky subtracted using the median of 
the dithered exposures, and then "shifted-and-added" using IRAF. They were also corrected 
for distortion 13 !. 

We achieved high-precision relative astrometry by geometrically transforming the pre- 
explosion images to match the post-explosion ones. We first "drizzled"0 the pre-explosion 
images for each band to the higher NIRC2+AO resolution. Using 5-7 point-like sources 
in common between the two datasets and the IRAF tasks geomap and geotran, we carried 
out a geometrical transformation between the two sets of images. The positions (and their 



12 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the AURA, Inc., under cooperative agreement with the National Science 
Foundation (NSF). 

13 http:/ /www2. keck. hawaii.edu/inst/nirc2/forReDoc/post_observing/dewarp/. 
14 http:/ /www. stsci.edu/hst/wfpc2/analysis/drizzle. html. 
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uncertainties) of the SN and the progenitor candidate are derived by averaging the mea- 
surements from two centroiding methods, the task daofind within IRAF/DAOPHOT and 
imexamine within IRAF. The differences between the SN and the progenitor candidate po- 
sitions, compared with the total estimated astrometric uncertainty, are given in Table [TJ As 
an additional check, we transformed the NIRC2+AO image relative to only the WF3 chip 
image in F814W at its native resolution (O'.'l pixel -1 ) and performed the registration. In this 
case, the positional difference is 19 mas, within ~ la of the uncertainty in the astrometric 
solution. Note that no other source was located within a 5a radius from the progenitor 
candidate position identified by Li et al. (2009); the Hll region to the northeast is located 
at > 25a (see Fig. ED. 



4. The Nature of the Progenitor 



We also measured photometry of the HST images using HSTphoO (IDolphin 1 120001 ). 
The output from this package automatically includes the transformation from F555W and 
F814W to V and I. 

Adopting a distance modulus to NGC 1832 derived from the recession velocity corrected 
for Virgo infall (32.09 ±0.30 mag; this uncertainty arises from a possible 250 km s _1 peculiar 
velocity 16 !) and the assumed SN extinction, we find that the absolute magnitudes of the 
progenitor candidate are M v = —7.80 ±0.33 and M® = —8.75 ±0.32, entirely consistent with 
a highly luminous supergiant. The intrinsic color, (V — I)q = 0.95 ±0.21 mag, is significantly 
more "yellow" than the colors of normal red supergiants (RSGs; e.g., iDrilling fc Landolt 
2000l ). 



We can directly determine the metallicity in the SN environment from the CFHT high- 
resol ution spectrum, by m easuring; the Ha and [N II] A6584 line intensities and applying 
the (jPettini &: Pagel 1 120041 ) cubic-fit relation between the ratio of these two lines and the 
oxygen abundance; no correction is made for the low extinction. Doing so, we find tha t 12 



+ log(0/H) = 8.67. Given that the solar value is 8.69 ± 0.05 flAsplund et al. Il2009h . we 



consider it most likely that this environment has roughly solar metallicity. 

The star's color corresponds to an effective temperature T eS = 5300 ± 500 K and a V 
bolometric correction in the range —0.47 to —0.13 mag, for an assumed surface gravity log g = 



lD HSTphot is a stellar photometry package designed for use with WFPC2 images. We used vl.l.7b, 
updated 2008 July 19. 

16 From NED, http://nedwww.ipac.caltech.edu/ 
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+0.5 (Kurucz Atlas 9 models, CD-ROMs 13, 18). This results in L bo i = io( 512±0 - 15 ) L 
(assuming M bo i(©)= 4.74 mag). In Figure [3] we show a Hertzsprung-Russell (HR) diagram 
i ncluding the progen itor candidate. In addition, we illustrate model evolutionary tracks 
(IHirschi et al. 1 120041 ) for stars with M ini = 15, 20, and 25 M , with rotation (t> ini = 300 



kms 1 ) and without rotation. 

The location of the progenitor candidate in the HR diagram is clearly not consistent 
with the 15 tracks, an initial mass which lies within the range for SN II-P progeni- 
tors (ISmartt et al. II2009I ). Unfortunately, the mass bins for these tracks are large (5 M 
increments). However, taking into account the uncertainty in the progenitor candidate's 
luminosity, interpolating by eye between the tracks implies that M- m i « 18-24 M . We 
note that this is con sistent with the upper limit of < 20 M on the SN 1980K progenitor 
(ISmartt et al. 1120091 . using more curre nt theoretical tracks) and the lower limit of ^ 17-18 
M on the SN II-L 1979C progenitor (jVan Dvk et al. Ill999h . 



That such a massive progenitor would be somewhat bluer than would be expected for 
the normal RSG progenitors of SNe II-P is also consistent with the star being in a post-RSG 
phase. Such an expectation follows from the theoretical models by Hirschi et al. (2004) 
for rotating stars at higher masses; the rotating models are more luminous and evolve to 
the RSG phase before the ignition of He burning, which results in higher mass-loss rates, 
the loss of most of the H envelope before the termination of He burning, and evolution 
toward the blue before the terminal points. From Figure [3] one can see that rotation makes 
little difference in the evolution of a 15 M star; both the rotating and nonrotating tracks 
terminate as a RSG. However, for a 20 M star, rotation has a profound evolutionary effect. 
Though we cannot tell from these models what occurs for stars in the 15-20 M range, it is 
reasonable to assume that rotation begins to affect those models nearing 20 M . Thus, we 
might expect the progenitor candidate, based on our mass estimate above, to have evolved 
toward the blue before explosion. 

Although more luminous and hotter than the SN 2009kr progenitor candidate, similar 
behavior has occurred for the post-RS G star IRC+10420, a m ass-losing hypergiant that is 
transiting the so-called "yellow void" ([Humphreys et al.ll2002l ). Likely more analogous are 
the "anomalous" yellow supergiants in the Magellanic C louds, which are in the mass range 
15-20 M and show evidence for post-RSG evolution (IHumphreys et al.lll99ll ). albeit at 
subsolar metallicity. 



Of course, there are caveats; as we found for the peculiar SN II-P 2008cn (jElias-Rosa et al 



20091 ). a yellow progenitor color could result from evolution in an interacting binary. Ad- 
ditionally, although the difference between the SN and progenitor candidate positions are 
practically within the total uncertainties, these are still large, and, given the host galaxy's 
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distance (a single NIRC2 pixel corresponds to ~ 5 pc), it is possible that we have not 
identified the SN progenitor at all, but rather a compact star cluster. 

If so, adopting the SN extinction, we can fit the V and / fluxes for the progenitor can- 
didate with the Starburst99 cod model spectral energy distributions and find an excellent 
fit with a 10 Myr track (Fig. [4]). If the SN progenitor were a member of this cluster, this 
age estimate would still be consistent with the lifespan of a ~20 M star. The 8 Myr model 
also provides a good fit. However, the lowest possible cluster luminosity that exceeds the 
luminosity of the most massive 8 Myr-old s tar is My « —8.1 mag, and for 10 Myr it is 
My m —7.6 mag (jCervino &: Luridianal 12004 ). Both of these are roughly consistent with the 
candidate's luminosity. Furthermore, we would ordinarily expect compact star clusters to 
have My < —8.6 mag (IBastian et al. 1120051 ; ICrockett et al. 1120081 ). Additional dust obscu- 
ration, beyond our assumed extinction, would imply that a putative cluster could be more 
luminous, but it would also be bluer and younger. 



A similar study has been done contem poraneously by Fraser et al. I (120101 ). identifying 
the same candidate progenitor. However, iFraser et al. I ( 120101 ) suggest that the SN is a 
spectrally peculiar SN II- P, with only r-band photometry reported. Differences also exist 
in the metallicity estimate and theoretical models employed, leading to somewhat different 
estimates for M ini (15+4 M o bv lFraser et al. Iboioh . 



5. Conclusions 

Based on ~ 3 months of follow-up observations from discovery, we conclude that 
SN 2009kr is a SN II-L, a relatively rare subclass of CC SNe. We identified an object in pre- 
explosion HST images which agrees astrometrically with the SN 2009kr position, using NIR 
SN images obtained with NIRC2+AO on Keck. From the SN distance and extinction we find 
that the object is consistent with a highly luminous supergiant star. Placing the inferred 
Lboi and T eS of the progenitor candidate on an HR diagram, we infer M ini « 18-24 M . Its 
yellow color implies that the star may have exploded in a post-RSG evolutionary state, which 
is predicted by our assumed theoretical stellar models. If this star is the SN 2009kr pro- 
genitor, this would be the first-ever direct identification of a SN II-L progeni tor. The star's 



mass estimate also may be a link between the upper mass ra nge for SNe II- P (ISmartt et al. 



20091 ) and the estimated progenitor mass for SN Iln 2005gl (jGal-Yam fc Leonard 1120091 ). 



Ultimately, the definitive indication that we have identified the SN 2009kr progenitor is 



17 http://www. stsci.edu/science/starburst99/, Vazquez fc Leitherer < 20051 ) 
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through very late-time imaging, such as wit h HST (IMaund fe Smarttll2009i ). However, as a 
SN II-L, akin to the long-lasting SN 1979C ( Milisavljevic et al. 20091 ). we may have to wait 



at least a decade to image the SN site using the James Webb Space Telescope. 
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SN 1992H (II — P) 
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rest wavelength (A) 

Fig. 1. — (a) SN 2009kr optical and NIR light curves, together with those of other SNe II-L 
and SNe II-P. The comparison light curves are adjusted in time and magnitude to match 
those of SN 2009kr. (b) SN 2009kr spectral sequence, along with comparison spectra of 
other SNe II-L and SNe II-P. All spectra have been corrected for their host-galaxy recession 
velocities and for reddening (see text). Ages are relative to B maximum light. The continuum 
of the t ~ 116 day spectrum of SN 2009kr appears to be significantly contaminated by light 
from the nearby star cluster. 
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Fig. 2.— (a) Subsections of the pre-explosion HST+ WFPC2/F814W images of NGC 1832, 
and (b) the post-explosion AO K p image of SN 2009kr with Keck II+NIRC2. The approx- 
imate positions of the candidate progenitor and the SN are indicated by tick marks, (c) 
5" x 5" detail of the HST image. The 3a positional uncertainty ellipse is 0'/03 x 0'.'09 in 
radius (see Tabled]). 
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Fig. 3. — An HR diagram showing the Lboi and T e g for the candidate progenitor of SN 2009kr 
(filled circle). Model stellar evolutionary tracks for solar metallicity (Hirschi et al. 2004) are 
also shown for a rotation of v ini = km s _1 (dot- dashed, solid, and dashed red lines) and 
fini = 300 kms" 1 (dot-dashed, solid, and dashed blue lines). 
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SN 2009kr + Starburst99 best-fit models 
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Fig. 4. — Comparison of the progenitor candid ate (filled circles) with mo del spectral en- 
ergy distributions for star clusters (Starburst99; I Vazquez fc Leitherer 112005c lines) for solar 
metallicity. The model parameters (age and total mass) are indicated in the legend. 
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Table 1. SN and Progenitor Candidate Position Comparison 





F555W (a/5) 


F814W (a/5) 


Uncertainty in the progenitor position (mas) 


0/2 


0/24 


Uncertainty in the SN position (mas) 


8/13 


8/13 


Geometric transformation (mas) 


4/15 


5/13 


Total uncertainty (mas) 


9/20 


9/30 


Difference in position (mas) 


13/4 


7/22 



Note. - - Uncertainties (la) in the SN and progenitor candidate right ascen- 
sion (a) and declination (5) positions for each band, in milliarcsec (mas), were 
estimated as the standard deviation of the average measurements. Uncertainties 
in the geometric transformation were derived from the differences in the fiducial 
star positions before and after the transformation. The total uncertainty is the 
quadrature sum of all uncertainties. The last row lists the residual difference 
between the SN and progenitor position after the geometric transformation. See 
text. 



